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We study the effect of the cr(600) and ai(1260) resonances in the p° — » 7r + 7r~7 decay, within 
the meson dominance model. Major effects are driven by the mass and width parameters of the 
<r(600), and the usually neglected contribution of the ai(1260), although small by itself, may become 
sizable through its interference with pion bremsstrahlung, and the proper relative sign can favor the 
central value of the experimental branching ratio. We present a procedure, using the gauge invariant 
structure of the resonant amplitudes, to kinematically enhance the resonant effects in the angular 
and energy distribution of the photon. We also elaborate on the coupling constants involved. 
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I. INTRODUCTION 

The p° meson can be produced via e + e~ annihilation. 
Dedicated experiments like KLOE, SND and CMD-II 
have performed extensive studies on the decay modes of 
this meson 0, 0, H, The measured branching ratio 
of the p° — > 7r + 7r~7 radiative decay is 9.9 ± 1.6 x 10~ 3 
The theoretical approach is able to explain the ex- 
perimental result, at the current precision, considering 
the radiation off pions alone. This can be estimated in 
a model independent way, predicting a branching ratio 
around one standard deviation above the central value 
[fil> B 0- The maximum energy that the photon can 
carry out in this process is of u> = 334.8 MeV, which 
is large enough to suggest that contributions beyond the 
soft photon approximation may be relevant in future ac- 
curate measurements. 

The introduction of the subleading contributions is 
model dependent. If the scalar contribution is identified 
with the cr(600) resonance then both mass and width may 
be obtained by making a fit to the data. Still, the range 
of values obtained using this procedure are not as pre- 
cise as the determination from other methods [1, |j| . The 
other allowed contribution for the decay is of the axial- 
vector form. The coupling of the relevant amplitude is 
estimated in the chiral perturbation theory to be propor- 
tional to the parameters L\ + L r w ~ 1.4 x 10 -3 [10| and 
therefore is commonly disregarded in the analysis 7] . On 
the other hand, meson dominance identifies this contri- 
bution with the ai(1260) state which, at first, is far from 
being on-shell and then is neither taken into account. 

The improvement on the experimental precision will 
make it possible to distinguish the resonant contribu- 
tions, and therefore a clear theoretical estimate will be 
important. In particular, we can ask whether such con- 
tributions can help to bring the theoretical predictions 
closer to the central experimental value and if there is a 
window where the resonant contributions are more rel- 
evant. The quark structure and proper characteriza- 
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FIG. 1: Feynman diagrams that contribute to the decay 
p° — > 7r + 7r~7. 



tion of such states is also a matter of ongoing research 

In this article we use the meson dominance model to 
provide a full account of these resonances. We perform 
a self-consistent determination of the effective coupling 
constants involved by reproducing observables in other 
decays, without invoking to one-loop corrections. Since 
the radiation is dominated by pion bremsstrahlung, the 
resonant contribution upon its interference with it can 
become relevant. Here we compute the contributions to 
the branching ratio from the different sources which may 
also provide a hint to the relative sign of the ai(1260) am- 
plitude by favoring the experimental central value. The 
di-pion invariant mass spectrum and the angular and en- 
ergy distribution of the photon are used to look for sensi- 
tivity to the resonances. We exploit kinematical configu- 
rations where for the latter may lead to enhancements of 
the resonant contributions. The cr(600) mass and width 
effects are explored throughout the analysis. 
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II. 



DECAY AMPLITUDES 



Let us state the conventions for the process p°(q, if) — » 
ir + (p )ir~ (p) 7(fc,e); q, p , p and fc are the correspond- 
ing 4-momenta and, 77 and e are the polarization tensors 
of the vector meson and the photon respectively. The 
Feynman diagrams contributing to the decay are shown 
in Figure (JTJ) . Diagrams and [TJ: are model in- 

dependent and correspond to the Low amplitude [Til ]. 
Diagrams QJi, [JJj and QJ are model dependent and, in 
the meson dominance model, include the intermediate a 
and cii resonances. We will use the dependence on these 
states to label the corresponding amplitudes. The inter- 
mediate state /(980) is not taken into account since it 
is considered to be dominated by the KK channel [l3| . 
Therefore, the total amplitude can be written as: 



Mr. 



Ml + M„ + M n + + M n - , 



(1) 



III. SQUARED AMPLITUDES AND 
INTERFERENCES 



We choose p ■ k and p ■ k as the independent variables 
to make the dependence on the photon energy explicit. 
We list below some of the relevant squared amplitudes 
and interferences: 



/ M 2 
16nag p ^L m ff - — 



M 2 p 



K - ml 



M 2 p p ■ k 
p ■ k m\ 



p ■ k 



(7) 



2p-k , (8) 



where the Low amplitude [f| Ml = M a + Mb + M c , 
can be split into two parts, owing to its dependence on 
the photon energy: Ml = Mei^ 1 ) + A^o( w °), where 



M e = 2ie g p7T7r p ■ r) L ■ c* 



(2) 



2Re M e M a = 16na g p ^ g pai g av7r p ■ kp ■ kL 2 x (9) 
/ M 2 p -m 2 a -2q-k \ 

I [M 2 p -m 2 a -2q- k} 2 + m%V% ) ' 



M = 2ie g p7! 



k ■ r\ 



p ■ k 



(3) 



and V = p v jp- k—p' v ' jp' ■ k satisfies L- k = 0. Therefore 
Ml is explicitly gauge invariant. g p7r7r is the coupling 
constant between the rho and the pions, as indicated by 
the subindex. The amplitudes for figures [TJi, and [TJ 
are gauge invariant by themselves and are given by: 
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tcgpcrj go 
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(p + p ) 2 - ml + im^Ta 



, (4) 



2Re M e M 



= 4e g plT1T g ai7r7 g aiP ^P ■ kp ■ kL 2 x (10) 
[ml + p ■ k] [ml - m 2 al + 2p ■ k]\ 



[ml + 2p ■ k - m 2 al } 2 + m 2 al T 2 al ' ' 



2Re M e M^ _ = 4e g plxlx .g QlTr7 g aiP ^V • kp ■ kL 2 x (11) 
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(6) 



respectively, where g a -Kn, gat^-y and g ai piv are effective 
coupling constants. M a (T a ) and m ai (T ai ) are the 
corresponding particle masses (full widths). The reso- 
nant propagators have been assumed to be of the com- 
plex mass form 15j, which corresponds to replacing 



m 2 — > m 2 — imT in the non-resonant propagator. 



We have not written the square of the model depen- 
dent amplitudes and the remaining interferences but they 
are actually taken into account in the calculation. The 
Low interferences of order w _1 are null in accordance 
with the Burnett-Kroll theorem [Hl | and a term propor- 
tional to L 2 from the M e M.\ interference was absorbed 
into eqn. 0. By inspection of eqns. ((71 - [TTj) we ob- 
serve that, in addition to |A"f e | 2 ) all the interferences are 
proportional to L 2 , a property we showed in a previous 
work to hold whenever there is an interference between 
the electric charge radiation and any gauge invariant am- 
plitude [ijj • This property will allow us to kinematically 
enhance the model dependent contribution by properly 
choosing the region where L 2 is maximum |18( . Although 
this enhancement is also promoted to the dominant elec- 
tric charge contribution, the latter receives a natural sup- 
pression as the photon energy increases, owing to the ui~ 2 
dependence, while the model dependent ones are of order 
u° and higher. 
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M a 


r CT 


Sigma 


Low + a 


Tot. (+) 


Tot. (-) 


MeV 


MeV 


10 -5 


10~ 3 


10~ 3 


10 -3 


500 


500 


1.7 


11.56 


11.64 ±0.03 


11.49 ± 0.01 


500 


450 


2.0 


11.57 


11.64 


11.50 


500 


400 


2.4 


11.57 


11.64 


11.50 


500 


350 


2.9 


11.58 


11.65 


11.50 


450 


500 


-5.8 


11.49 


11.56 


11.42 


450 


450 


-6.1 


11.49 


11.56 


11.42 


450 


400 


-6.4 


11.48 


11.55 


11.41 


450 


350 


-6.5 


11.48 


11.50 


11.41 


400 


500 


-11.7 


11.43 


11.49 


11.36 


400 


450 


-12.4 


11.42 


11.49 


11.35 


400 


400 


-13.0 


11.42 


11.49 


11.35 


400 


350 


-13.7 


11.41 


11.48 


11.34 


350 


500 


-17.0 


11.38 


11.45 


11.31 


350 


450 


-17.9 


11.37 


11.44 


11.30 


350 


400 


-18.8 


11.36 


11.43 


11.29 


350 


350 


-19.7 


11.35 


11.42 ±0.05 


11.28 ± 0.07 



TABLE I: Branching ratios from several contributions for a 
set of values of the a parameters and a cut off on the photon 
energy of 50 MeV . Low = 11.547 x 10~ 3 . To have an idea 
of the effect from the uncertainties in the coupling constants 
and ai mass and width, we have included error bars in the 
first and last row. See text for details. 



IV. RESULTS 

In order to make an estimate of the observables, we 
first address the problem of finding the proper values of 
the coupling constants. The g plT7T and couplings can 
be written in terms of masses and widths as: 
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FIG. 2: Di-pion invariant mass due to the pions 
Bremsstrahlung for the p° — > 7r + 7r _ 7 decay. 



which lies in the region of |6 to 7\, and produces an over- 
estimate of the branching ratio for p° — > an — > 7r°7r°7 
[2lj . Our ecmation UU is valid in the limit of T a = 0. In 
reference [22j the coupling was computed including the 
width effect and relying on the experimental value for 
T(p° — > o~tt — > 7r°7r 7 ) = {2.9±H ± 0.6) keV [|. Using 
theirs master equation and varying m a and T a from 350 
to 500 MeV produces g pai = 0.22 to 0.89, which is of 
the same order than our value, although in a different 
approximation. Therefore the use of either values will 
produce similar results in the range we are exploring. 
An estimate of the g aip7T coupling can be obtained from 
the a\(q,r]) — > p(k,e)ir decay. The amplitude in its sim- 
plest on-shell form becomes [231 ] : 

q ■ e* 

M(a! -* pn) = fa lP ^(v ■ e* - k ■ -q— ) (15) 

Our coupling is related to this by g aip7T = —faipir/k ■ q. 
Assuming that the full width is dominated by the pir 
channel, the coupling is 



f 



127TT<r oi 



[p-7r°; 

where [pit] = [1 ± 



[p°7T- 



(3.7 to 5.9 GeVf (16) 



where we have used for definiteness T p = 150.7 MeV and and k ■ q = {m 2 ai ± Mj - m£)/2, and we have used 

m 2 ai = 1230 ± 40 and T ai = 250 - QOOMeV. This 
compares well with the prediction of f aip7r — 4.8 GeV in 
the quark model [231 ] . 

Finally, using vector meson dominance ar- 



m p = 77hMeV. 

For the g pcr ~ ( coupling, we use the expression depending 
on the radiative width T ptTJ whose value was estimated 
to be r p(T7 = 0.23 ± 0.47 keV or 17 ± 4 keV [T3]. 



2 _ 3 
Spry n L P 



2M n 



Ml 



(14) 



guments ;24| , we 
g ai TT-y as follows: 
7p = 2ay/KM p /3r(j>P 



can 

gaiTZ~f 



relate the 



and 



la^hp, where 



o+e-) = 5.012. 



Taking the prediction for r pCT7 = 17 ±4 keV and varying 
m a = 350 - 500MeV produces g 2 ptn = (1.82 to 7.46) x 
10" 7 MeV- 2 . To compare with the one obtained in 
[20T ] we re-parameterize the coupling by a M p factor ( 
5po-7 —* fjpa-y = Mpgpa-y — > 0.33 to 0.67). Therefore, it is 
one order of magnitude smaller than the extracted in (20| 



Errors in the observables coming from the uncertainties 
on the coupling values will be specified below. 

The radiative decay width can be computed using the 
expressions for the squared amplitudes and interferences 
separately. In table HI we present the different contri- 
butions to the branching ratio for a set of values for 
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the mass and width of the ct(600). We have introduced 
a 50 MeV cut in the photon energy to avoid the in- 
frared divergence, then its allowed region is: uj cut < 
u < (Mp ~ /2M P = 334.8 MeV. The upper value 
requires to go beyond the soft photon approximation, 
given by the Low amplitude. It is worth to mention 
that the higher the cut in the photon energy the bet- 
ter the sensitivity to the resonant parameters. Here we 
stick to this value to compare with the available data. 
The Low amplitude contribution to the branching ratio 
is Low = 11.547 x 10~ 3 , about one standard deviation 
above the experimental value of 9.9 ± 1.6 x 10~ 3 [l|, for 
the same cut off. The column labeled "Sigma" corre- 
sponds to the contribution from the a amplitude itself 
and its interference with the Low amplitude. This be- 
comes larger for smaller values of the a parameters and 
can even flip the sign. Column labeled "Low +sigma" is 
the sum of the previous column plus the Low contribu- 
tion. Columns labeled "Tot(+)" and "Tot(-)" correspond 
to the total branching ratio when the a\ interferences 
with the Low amplitude is also included. Depending on 
the sign of the a\ amplitudes the contribution can be ei- 
ther oi(+) = 6.8 x 10" 5 or oti(-) = -6.6 x 1CT 5 . 

From table HI we observe that the inclusion of the reso- 
nances can help to bring the prediction closer to the cen- 
tral experimental value. In particular the sign on the a\ 
amplitude can either improve or worsen the agreement. 
Still, the major effect is driven only by the m a and r CT 
parameters. To have an idea of the effect from the un- 
certainties in the coupling constants discussed previously 
and Oi mass and width, we have included error bars in 
the first and last row of Table Q] The contributions from 
the cti itself is of order of 10 -6 and neglected in the re- 
sults. 

The di-pion invariant mass is interesting on its own 
and the experiments usually report it as the main ob- 
servable. However, in our case the individual effects from 
the interferences between the electric radiation and the 
resonances are very mild, as can be expected from the 
results for the branching ratios. Just for illustration, in 
Figure[5]we plot the corresponding Low amplitude, which 
dominates this observable. 

The fact that all the leading interferences are propor- 
tional to L 2 , allows to look for an enhancement of the 
effects from the resonances by choosing the kinematical 
configuration where L 2 is maximum. In the p° restframe 
it can be written as L 2 oc 1 — cos 2 8 where is the angle 
between the photon and n~ 3-momenta. In addition, the 
dependence on the photon energy of the interferences is 
of order uj° while the dominant Low contribution is of 



order 



This suggests that an appropriate observ- 



able could be the angular and energy distribution of the 
photon. 



dT = £ 

E=E+,E- 

?± 



\M\ 2 M p 



(27T) 



F\E) 
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FIG. 3: Photon energy spectrum due to the Low emission 
for 8 = 85° and 5° (solid and dashed lines respectively). 




where, E ± labels the roots of F(E) 



M 2 - 2M p E 



0.1 0.15 0.2 0.25 0.3 0.35 
Photon energy, x=w/M p 

FIG. 4: Photon energy spectrum due to the interference 
between electric and a emission for 6 = 85° and 5° (solid and 
dashed lines respectively). Here m CT = r CT = 400 MeV . 



2M p u> + 2(Ew — pujcosO), in the p° restframe, and we 
have introduced the dimensionless variables y = cosQ 
and x = uj/Mp. The maximum value for uj is (M 2 — 
2m^Mp)/2(M p - In Figures H H and [5] we plot 

dT/T nr dydx (normalized to the non-radiative width, 
r nr ) for two angles of the photon emission, = 85° 
and 5° (solid and dashed lines respectively). Figure 
[3] includes the Low contribution, Figure [J] the inter- 
ference between Low and cr, and Figure [5] the inter- 
ference between Low and a\ emission. Here we have 
used M a = r CT = 400 MeV. We can observe that, al- 
though small, the contributions from the resonances can 
be strongly enhanced by choosing the proper angle off 
emission, about 50% for the a and upto 85% for the a\ 
in the current set up of relative angles. This enhance- 
ment is not promoted to the Low emission at the same 
proportion, which is mildly affected and even suppressed 
for large values of the photon energy as expected (Fig- 
ure [3]). The total photon energy and angular spectrum is 
certainly dominated by the Low emission but this is free 
of relevant theoretical uncertainties and can be safely re- 
moved from data. 
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FIG. 5: Photon energy spectrum due to the interference 
between electric and ai emission for 6 — 85° and 5° (solid 
and dashed lines respectively). 

V. CONCLUSIONS 

We have studied the p° — ► n + ir~-f decay in a 
self-consistent approach based on the meson dominance 
model, and included both the cr(600) and ai(1260) reso- 
nances. We determined the corresponding coupling con- 
stants involved and, in particular, the value for g pcr -y is 
one order of magnitude smaller than the one extracted in 
previous studies [2C| but similar to the estimation made 
in [22l | for the particular range of parameters we did ex- 
plore. Rigorously one must use the coupling from the 
latter where T a is taken into account. 

We identified the different contributions to the branch- 
ing ratio, dominated by the pion bremsstrahlung whose 
contribution alone lies one standard deviation above the 
experimental central value (measured with a precision 



of about 16%). The resonant contributions upon in- 
terference with the Low radiation, although small, can 
be of relevance for future accurate measurements of the 
branching ratio and be sensitive to the resonance param- 
eters. In fact, this would provide a hint on the relative 
sign of the axial amplitudes by requiring the theoretical 
prediction to lie closer to the central experimental value. 
In order to distinguish the effects coming from the <r(600) 
parameters a precision smaller than 5% is required, while 
to be sensitive to the a\ parameters at least the one per- 
cent level is required. 

By another hand, the di-pion invariant mass spectrum 
was computed and shown to be saturated by the pion 
Bremsstrahlung . 

Exploiting the structure of the leading interferences, 
we tuned a kinematical configuration where resonant con- 
tributions can be enhanced, namely the photon angu- 
lar and energy spectrum. In particular the effects can 
be enhanced for quasi-transversal emission compared to 
quasi-collinear emission of photons, with respect to the 
7T _ 3-momentum. Our treatment is useful for looking 
for enhancements of the resonances in decays of the form 
p° — » 7r + 7r~7, since it exploits the radiation structure of 
the external charged particles and can serve as a comple- 
ment to estimates from decays of the form p° — > 7r°7r°7 
which are mainly driven by model dependent contri- 
butions and where charged particles contributes only 
through loops and therefore our approach can not be ap- 
plied. 
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